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Abstract 

In many theories of reheating starting from the classical spatially ho- 
mogeneous inflaton field, its accompanying inhomogeneous part (which 
arises from primordial quantum fluctuations) is treated as a first or- 
der perturbation. We examine some consequences of treating it non- 
perturbatively in a model where a first order treatment is invalid. In 
particular we consider effects on the long-wavelength curvature pa- 
rameter £. 

1 Introduction 

Starting from minimal models of inflation with one scalar field, <p, the ad- 
dition of a scalar field, \i interacting with it has been much investigated, 
particularly in the context of parametric resonance |L|, 0, [|, |], §, |6|, [7], §. The 
inflaton field, <p, is equal to <po + y?i, where <po is the classical field giving rise 
to inflation through a potential V(<p) and <pi arises as a primordial quantum 
perturbation, (fx is associated with perturbation to the metric also arising 
from primordial quantum perturbations. In longitudinal gauge (and in the 
usual case of absence of space-space off-diagonal elements of the stress-energy 
tensor) the FRW metric appears with the perturbing field ip as 

ds 2 = a(r) 2 (l + 2ip)dr 2 - a(r) 2 (l - 2^)5 ij dx i dx j . (1) 

The equations of motion imply that tp\ and ip arise having identical k- 
component quantum creation and annihilation operators. In this minimal 
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model the small k (long wavelength) components of ipi and if> are responsible 
for the observed cosmic microwave background radiation fluctuations (CM- 
BRF). Though an additional field \ could also have an influential classical 
part, thus forming a two-field inflation model we shall consider single field 
inflation with \ only arising from its own primordial quantum fluctuations. 
In fact with the well-used interaction g 2 <p 2 x 2 0> @ ; [| 0; 3 there are reasons 
for any x field being strongly suppressed before reheating begins. So in this 
paper on reheating we have initially three quantum fields tpi,ip and x >the 
latter being composed of different quantum operator components to those of 
(pi and if), which are identical. (The x field fi as a seemingly negligable value 
at the beginning of reheating compared with the other two but can increase 
to logarithmically comparable values through resonance effects, as we shall 
see.) 

These quantum fields decohere into classical stochastic fields. There are a 
number of treatments of decoherence but the work of Polarski and Starobin- 
sky || , giving a gradual decoherence with the expansion of the universe and 
the multiplication in the number of particles, is particularly relevant. The 
stochastic fields inherit, in transmuted form, important properties of their 
quantum progenitors; in particular ensemble averaging of the products of 
stochastic fields gives similar results to taking vacuum expectation values of 
corresponding products of quantum fields. 

For example, arising from the interaction g 2 ip 2 x 2 , the equation of motion 
of the ip Q field contains a term g 2 ifioX 2 and the only simple way to deal 
with this is to take an ensemble average of the x 2 term, which we denote 
by < X 2 >• I n terms of the mode functions this gives exactly the same 
expression as the vacuum expection value, < x 2 >o, of the corresponding 
quantum fields. Thus, since the mode functions themselves are smooth, there 
is a smooth transition from the quantum era to the classical stochastic era. 

A very significant use we make of ensemble averaging is in the calculation 
of the Hubble parameter, H(a'/a = aH), The picture during the inflation 
period is that the classical inflaton field dominates the pressure and energy 
densities and gives H. However it can happen during reheating - and does in 
the examples we shall treat - that with the diminution of (fo the contribution 
of (fx to the energy density rapidly becomes larger than that of <po, and 
its contribution to the spatially homogeneous Einstein equations and the 
calculation of H cannot be ignored or treated as a perturbation. Thus if 
is the Einstein tensor evaluated from the unperturbed (ip = 0) metric, 
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ds 2 = a{r) 2 dr 2 — a{r) 2 8ijdx l dx\ we write 



.Gf =< T? > . (2) 

As one example this means that a term (<p' ) 2 appearing on the RHS of Eq.(0) 
receives an addition < (f[) 2 > + < (x'i) 2 >■ Thus the homogeneous energy 
and pressure densities, po and po, during reheating and consequently the 
magnitude of H, may be largely governed by such additions. 

As an illustration of some results of this approach in our model we eval- 
uate the well-known parameter ( [|K], |Tl| which, as a function of the wave 
number k, is defined [F2 through the /c-component, ip k ,of the metric pertur- 
bation ip 

( k = ^{H~ 1 ^ k + ^ k )/{l + w)+^ k (3) 

For those small values of k corresponding to the observed CMBRF we have 
that k 2 / a 2 -C H 2 for any H in the reheating period and thus ( k ~ — lZ k where 
lZ k is the curvature perturbation [13]. For such k and with the assumption of 
adiabaticity ( is constant through reheat and to the matter era [10, 111], [12L 



[TSJ, [TJ], [T^, So it has often been used as a tool in calculating the CMBRF 
in various models of inflation; there has been considerable interest in the 
effect of parametric resonance on this issue [16|, [17|, [Tj|. We shall examine 
the variation of ( both analytically and numerically in the context outlined 
above. 



2 Decay of the Inflaton and Preheating 

In this paper we shall not treat true reheating (otherwise defrosting) in which 
there is conversion in large part to a thermal (relativistic) fluid. Rather we 
shall deal with a first stage of the conversion of the classical field ipo into 
particles of the fields x an d <fi- This, specially with respect to Xi has been 
named preheating and we shall see that the particles tpi can stand equally 
importantly in this respect. This being understood we shall continue to use 
the word reheating to include also this preheating stage. 

2.1 The equations of motion 

The equations of motion divide into two classes, the spatially homogeneous 
and the spatially non-homogeneous. The former are those for those two 
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variables which are functions of time only, a(r) and v 9 o( r )- The latter are 
those in the scalar fields <^i(x, r), x(x, r),t/){x, r) which are functions of space 
and time; these equations we shall express in terms of the Fourier component 
mode functions <^fc(r), X&(t), ipk(T~) in a way we shall describe below. 

We emphasize that this is not a division into non-perturbed and pertur- 
bation equations. The equations are to maximum order in Lp\ and x as we ll 
as ifQ. We do treat ip as a perturbation and we examine the validity of this 
approximation. 

We now specify the model we use more precisely [19[ . 

The Lagrangian is 

L = f d 4 x^{l^, a + lx a X,a - V(ip) - V( X ) - V int (ip, x)} (4) 



where ip = (po + <p\ and the longitudinal gauge metric is given by Eq. ([[]). 
This leads to the energy-momentum tensor 

T* = + X lfl X,u - + \x a X,a - V(<p) - V( X ) - V int }5» (5) 

We take the field potentials during reheating to be 

v{v) = \m\ 2 - v( x ) = \m\ 2 - v mt (p, x ) = \gVx 2 - (6) 

The potential V(ip) = \m 2 ip 2 is that used in chaotic inflation theory but is 
applied in this paper only in the reheating period; we use a different, but 
smoothly joining, potential in the inflationary era fl9|j . 
The classical stochastic fields are of the form 



^i( x > r )= / 7— r r[e(k)(^ k (r)exp(2k.x) + e*(k)<^(r) exp(-zk.x)], (7) 



— - T [d(k) Xk (r)exp(ik.x) + d*(k) x *(r)exp(-ik.x)] ) (8) 

(27TJ 2 



where e(k) and d(k) are time-independent separately 5-correlated Gaus- 
sian variables || such that, where (....) denotes the average, 

(e(k)e*(k')> = <d(kK(k')> = ^ 3 (k - k'); (9) 



with all other two- variable averages being zero. <pu( T ) an d Xk( r ) are the mode 
functions in which we express the equations of motion. The field ^(x, r ) is 
expressed as Eq.(|7|) with the same gaussian variables but the mode functions 
being ^ k (r). 

The properties of the gaussian variables given above result, for example, 
in the ensemble average 

(tp^x, r) 2 ) = (2tt)- 3 J <?k^i (10) 

(The originating quantum fields are expressed as Eqs.(^,|]) except that 
the gaussian variables are replaced by appropriate creation and annihilation 
operators; the vacuum expection values < .... > of products of quantum 
fields are unchanged in form from the ensemble averages < .... >; so ,for 
example, (<fi(x, t) 2 ) repeats the same form as Eq.flTtjp.) 

Integrals such as that in Eq. (|T0|) occur throughout the equations of motion 
and must be evaluated numerically Qand we have to adopt a finite range of 
wave number, k. If the integrals diverge as k — > oo then the upper limit of 
the k integration forms a cut-off which is the crudest way of dealing with such 
ultra-violet divergences. Two points of view may be taken on this. Firstly 
this may be considered equivalent to a renormalization procedure in mass 
and other quantities. Secondly the Lagrangian used may be considered as 
an effective Lagrangian which has absorbed extra degrees of freedom coming 
from supersymmetry which eliminate divergences at higher momentum, the 
cut-off representing this effect. In our numerical work we have generally taken 
a cut-off to correspond to a wavelength of H^ 1 so that k cuto ff ~ 2naH; this 
has been used in other work involving ensemble averages. We discuss further 
on the effect of different evaluations. 

We can now proceed to write down the equations of motion. As the above 
remarks explain these are the same in the quantum regime as in the classical 
stochastic regime. First we must consider Tjf. Defining the density and 
pressure homogeneous parts of the energy- momentum tensor as Ph(t) = (T®) 
and Ph(r)5*- = —(Tj) we find, after ensemble averaging, that 

Ph{r) = ^+a 2 m% 2 + <^>)+^ (11) 

Ph(r) = ^-a 2 m% 2 + (^ 2 >)-a^ 

(12) 

1 This is discussed in detail in ref. ||19|. 
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V = ^ + K> + (x' 2 ) 

- 2 2 i 2/ 2\ 

m = m + g (x /■ 
and we have the Friedmann equation 



(13) 
(14) 



{a'/af = (15) 

As in Eq. fllOD the averages are independent of x ensuring the spatially 
homogeneity of the R.H.S. of Eq. (|I~5|) . We also see from Eq.fllOD that Eq.([L5|) 
is equally valid in the quantum regime since the formalism ensures that the 
vev of the quantum operators equals the average over the stochastic variables. 

The other spatially homogeneous equation is 

ipl + 2(a'/a)(f/ + a 2 mVo - 4<VVi) ~ 4(^VVi) + 2a 2 rh 2 (^ 1 ) = (16) 

We note the combined reaction of the inhomogeneous fields tpi and ip on 
the homogeneous inflaton field through the ensemble averaging, similarly to 
that in Eqs. (0) and (0). 

There remain the 3 spatially non- homogeneous equations which we write 
in the /c-component form. These components are specified as the complex 
mode functions Xk, and ipk- Their wave number dependence, given by the 
succeeding equations, is only on k = |k| . 

Xl + 2(a'/a) X ' k + (k 2 + a 2 M 2 ) Xk = (17) 
where M 2 is a function of r given by 

M 2 (r) = M 2 + sVo(r) 2 + g 2 (<fi(x, r) 2 )o (18) 



(pi + 2(a'/a)(p' k + (k 2 + a 2 m 2 ) Vk - Aip' Q ijj' k + 2a 2 m 2 V ^ k = (19) 

Though not indicated these equations actually hold for each separate value of 
k; thus consistency of Eq.(pJ[) justifies our previous statements that ipk, the 
mode function of the metric perturbation, should be associated with the the 
same gaussian operators (or, in the quantum regime, with the same quantum 
operators) as ip^. 

/($ k 
— r[e(k)Vk(r) exp(zk.x) + e*(k)^(r) exp(-ik.x)] (20) 
(Z7T) 2 
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where since ^(x, r) is dimensionless ip^ has dimension (mass) - ?. This asso- 
ciates the metric perturbation with the inhomogeneous part of the inflaton 
field without assigning priority to either. But the stochastic variables of Xk 
are independent. Thus no terms in ip^ appear in Eq. (fl7|) ; they are forbidden 
through ensemble averaging.. 
The mode equation for ip is 

V>fc + W k (a'/a) + 4> k {2{a'/a)' + (a /a) 2 ) = AnGa 2 5p k (21) 

where 5pk(Sp k ) is the fc-component of the non- homogeneous part of the mo- 
mentum (energy): 



hk = \[-(v + (<Pi,i<Pi,i) ft + (x,*X,t)/3)V'fc + fWk ~ 2( Pk(i )( Pi) ~ a 2 m 2 po¥k}- 
er 

(22) 

In addition there is the time-space Einstein equation which acts as an equa- 
tion of constraint on the initial values. (We use it to fix the value of ip' at 
the beginning of reheating.) 

ij/ k + (a'/a)il>k = 4iiGip' Q ip k (23) 
2.2 Initial conditions 

We need to have the values of H, the fields and their time-derivatives at the 
end ofinflation to supply the initial conditions for the reheating equations 
of motion of the previous section, 2.1. We have chosen to use a specific 
inflationary model having the advantage that the solutions are analytically 
expressible. This is power-law inflation with an exponential potential: V = 
Uexp(—\(p) where U, A are constants. This potential of the inflation era 
stands in place of the potential V((p) = \m 2 Lp 2 of the reheating era but 
otherwise the Lagrangians are the same and in particular both have the 
interaction potential Vi n t(f, x) = \9 2{ P 2 X 2 - We have ensured the correct 
continuity by imposing the Lichnerowicz conditions |20(|, as well as using the 



equations of motion and constraint appropriate to each era at the boundary 
03 1 . 

An essential feature in maintaining the analytic form of power-law infla- 
tion is that the x-Aeld be negligable through the relevant era. This is ensured 
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by the field interaction term in the Lagrangian. A simplified version of the 
mode equation, Eq. (|T7|) , in the inflationary era is 



(a Xk )" + [k + a 2 M 2 + g 2 aV ~ a"/a}(a Xk ) = (24) 

and ifl is greater than, or of the order of, m 2 Pl through most of the inflationary 
era. Thus the term in square brackets is large and positive resulting in a 
quasi-periodic type solution for axu- The many efold increase of a during the 
inflationary era indicates an exceedingly small value for Xk at the beginning 
of reheat |21], |22"f , which may be as little as 10~ 50 of its initial value. This is 
an important qualitative feature of our initial conditions. We have taken the 

— 1 /2 

initial value of x to be \Q~ n m Planck where n is of the order of 30. (For such 
small mode functions in the beginning of reheat the transition to classical 
stochastic functions cannot yet be made and the quantum complex formalism 
should be retained.) 

For power-law inflation the scale factor a oc (tj — r) p and we are free to 
choose p by specifying A in the potential |19[ . In the numerical results quoted 



we have chosen p — —1.1 (p 1 corresponds to exponential inflation). 

2.3 The metric perturbation,?/; 

As noted above ip is treated perturbatively in the equations of motion and 
we check the validity of this, in each particular case used, in the following 
sense. We shall discuss ip as having developed into a classical stochastic 
field, Eq.(PD|). The indeterminancy represented by the variables e forces 
us to consider the average over the product of these variables so that we 
evaluate < (if>(x, r) 2 > the result being (ip(x,T) 2 ) = (27r)~ 3 / d?kip k {T)ip k (r)* 

the same for every value of x. We require that J < (ip(x, r) 2 > be small 
compared with unity, since the relevant metric coefficient is a(r) 2 (l + 2ip). 
Our viewpoint is that in any particular case satisfaction of this requirement 
forms sufficient justification for the perturbative approach, because the only 
way we can mount a comparison of the revised metric coefficient with a(r) 2 is 
when we consider the basic equations to be those in configuration space, and 
then indeed the revision is just a perturbation through all space-time[n|. In 



the cases shown in the figures y < T ) 2 > is less than 10 2 ' 5 throughout 
the range shown. 
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3 Variation of ( During Reheating 



Using the equations of motion of the preceeding section we can now analyse 
the behaviour of ( to see under what conditions it may be constant and what 
may cause it to vary. 

Multiplying Eq.(|[) by §-£^(1 + w ) an d differentiating with respect to cos- 
mic time we find after some manipulation and using Eqs.(|ll]) and (|12| ) 

Q J 

-H(l + w )( = 4 ! + HiP + 2H^j - $ + HiP) - \n(p h + p h ) (25) 



d_ } ( , \ _ VWo + (_Wh) + (xx) + (0i,»yi,«)/3 + (x,»X,i)/3 

dt [n[ph + Ph) ~ z <% + (vi) + (x 2 ) + (<Pim,i)/3 + <*,i*.i>/3 { ] 

where we have omitted the last terms in Eqs.(|ll]) and (12) as giving rise to 
terms of second order in ip. (Indeed including them would make no difference 
to the arguments we shall give.) 

We now compare the RHS of this equation with an equation for ip deduced 
from the Einstein equations. Besides Eq.(|21|) there is also the time-time 
equation in ip: 

V 2 V» fc - 3(ip' k + {a'/a)ifj k ) (a' /a) = AnGa 2 5p k . (27) 

Our concern is with ip k of wave numbers relevant to the CMBRF, so we can 
drop the first term in this equation. Then subtracting it from Eq.([H]) and 
using cosmic time we obtain 

^ + 7HtP + 2(H + 3H 2 )tfj = -47iG(2m 2 <p oVk + & k ) (28) 

£=A(K*>) + <4>) (29) 
Now eliminate (po by Eq . ([16|) and (p' Q by Eq.(^) to get 



4+(iJ-2^oM+8(^^i)/^o)^fc+(2ij-2^o/^o+8^(^^i)/^Q+47rGOV'fc = 

(30) 

Eq. (|30|) has some structure in common with Eq. ( p5|) . Thus putting all sec- 
ond order terms in the stochastic fields equal to zero both the RHS of Eq. (p5|) 
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and the LHS of Eg. (ftO]) reduce to the expression ip k + (H — 2cp /ip )ip k + (2H — 
2H<p /(p )ip k , implying that ( = 0. This essentially replicates the conditions 
under which Mukhanov et al.|T2j demonstrate the constancy of £. Under our 
extended equations it seems most unlikely that ( be constant except when 
ipo dominates and thus the quadratic terms, < ... >, can be neglected; in our 
model this is so at the beginning of reheat but then tp Q rapidly decreases and 
so we should expect ( to be no longer constant. Our numerical calculations 
do indeed show this expected behaviour of (. 

It is interesting that in our theory, where ipi is not treated perturbatively, 
then ( is still not apparently constant when we eliminate the x field so that 
X = in Eq . (p5|) and Eq . (|30|) , as is born out numerically. In this connection 
we note that one of the relevant ways of distinguishing entropic from non- 
entropic reheating is to evaluate — ^; for the reheating to be non-entropic 
this should be zero ||16|| . The expression we find in our formalism and with 
if i not treated perturbatively gives H — ^ 7^ even when x is n °t present 
in the theory. 

In the figures we have illustrated some examples from our numerical cal- 
culations. Our interest is in the qualitative features of reheat in a model 
compatible with known observations rather than in making precise compar- 
isons with data. So we are interested in values of the parameters that roughly 
supply the usual requirements of inflation and the magnitude of the CMBR 
fluctuations. The results we give are for the values m = 10~ 7 mpi, initial 
ifo = 0.3, p = —1.1, M/m = 0.02 and g/m = 2 x 10 4 . Throughout we quote 
dimensionful results and parameters in units such that H = c = G = 1. 

In Fig.l we show the energy densities of ipo and (fx; we see that in this 
model the energy density of ipi soon becomes comparable to that of <po and 
finally somewhat exceeds it; it is only at the beginning of preheating that 
if 1 can be treated as a perturbation of <^o- Fig- 2 shows how ( is initially 
constant, that is when (pi is much less than (fo, but then develops strong 
fluctuations and later increases. Fig. 3 shows ( developing resonance at 5 
efolds after the end of the inflation era but we note that (, as shown in Fig.2, 
has begun to fluctuate strongly while x is still negligable thus showing the 
influence of quadratic terms in (p. And in Fig.4 we give the ( that results 
when we have eliminated x from the equations, again showing the influence 



2 The figures in the previous versions of this paper were wrong; the numerical code for 
the solution of the equations had an elementary error. Some details of the coding used 



here are given in the authors' paper bcp-ph/0109218 
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of ipi and ip quadratic terms in the fluctuations of (. 



4 Summary and Discussion 

We have considered a single field (ip) model of inflation, with another scalar 
field (x) which is naturally quiescent during inflation and becomes active 
post inflation with the possibility of parametric resonance and preheating 
as the classical part, po, of the inflaton field decreases and oscillates. We 
adopted a method having stochastic variables naturally succeeding the quan- 
tum operators of the early inflation era and have developed equations for the 
post-inflation period by taking averages over the stochastic bilinear forms in 
the scalar fields. This enables a non-perturbative treatment of the stochastic 
scalar fields x an d p>i (where <p = ip + ipx) though the metric perturbation 
field, if), has to be treated perturbatively. 

In our results as <p decreases, the energy density of the p>\ field comes 
to somewhat exceed (by a factor of the order of 10) that of the po and also 
somewhat exceeds that of the resonating and much grown \ particle field. 
Even while this latter feature could be to an extent parameter dependent we 
consider a first order perturbative treatment of <pi to be inappropriate. 

Some features are illustrated by the behaviour of ( which is constant at 
first but then changes as the effect of the fields p\ and/or \ becomes no 
longer negligible in the other equations of motion. Our algebraic analysis 
showed that ( must change when we include quadratic terms in either ipi or 
X- Thus the necessary particle field (inhomogeneous field) concommitant of 
the classical inflaton field does by itself produce changes, independent of x- 

The question is not whether these sort of effects occur but rather how 
strong are they? 



The recent work of refs. [|1J| [16] finds that they are very small by express- 
ing the time variation of ( in terms of the non-adiabatic pressure variation. 
This comes using the reasonable physical assumption - based on the approx- 
imate validity of the Robert son- Walker metric - that after smoothing on a 
cosmological scale below the ones of interest the spatially inhomogeneous 
variables such as 5pk,Spk can be treated as perturbations. There is no evi- 
dent agreement between that expression for ( and Eq.(p5l) above which has 
no assumptions on 5pk,Spk being perturbative. Further differences arise be- 
cause of the use, in refs. |T3|, [L6| , of power spectra of x 2 an d C i n magnitude 



criteria; this gives rise to expressions such as Eq.(18) of |Tj|, different from 
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any in this paper, which are typified for example by Eq.(^) for 6pk- 

In this paper we have taken a high momentum ('ultra-violet') cut-off 
to make the loop integrals finite. In the cosmological context such ultra- 
violet cut-offs occur naturally as the non-zero lattice spacing in numeri- 



cal configuration-space calculations || [23|. They are also used in other 
momentum-space numerical calculations such as those of ref.[16] where, as 
in this paper, the cut-off used is H~ l . The value of H is roughly that of the 
preheating period. This value could be regarded as natural in this context 
even though a cut-off is the crudest, though still-used, way of regularizing in 
field theory. A marked change in its value (which acts in comparable degree 
both on the y?i integrals and on the x integrals) leads to different energy 
densities and effective masses, and to marked variations of the changes in ( 
and other quantities. 

Except for its consistent use of stochastic variables the theory and its 
parameters that we have used is neither artfully designed nor unusual. We 
consider that it points up the potential importance of the inflaton particle 
field in reheat and points to a degree of caution necessary in the choice of 
theory and parameters if that field is to be treated by first order perturbation. 

We thank Andrew Liddle, David Lyth and Luis Mendes for some conver- 
sations and Bruce Bassett for a question. 
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Figure 1: Logarithmic plot of energy densities of tpo (A)and tpi (B) versus 
the number of e-folds of expansion in the reheat era. In all figures densities 
are in units m 2 x m 2 planck , where m = 10~ 7 m P i anck . 
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Figure 2: Plot of lg(£) verss the number of e-folds of expansion in the reheat 
era. 



14 




4 5 
n(cfolds) 



Figure 3: Logarithmic plot of the energy density of x versus the number of 
e-folds of expansion in the reheat era. 
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Figure 4: Plot of lg(C) when the x field is switched off. 
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